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Abstract 

A detailed study is presented of the counterrotating model (CRM) for electrovacuum sta- 
tionary axially symmetric relativistic thin disks of infinite extension without radial stress, in 
the case when the eigenvalues of the energy-momentum tensor of the disk are real quantities, 
so that there is not heat flow. We find a general constraint over the counterrotating tangential 
velocities needed to cast the surface energy-momentum tensor of the disk as the superposition 
of two counterrotating charged dust fluids. We then show that, in some cases, this constraint 
can be satisfied if we take the two counterrotating tangential velocities as equal and opposite or 
by taking the two counterrotating streams as circulating along electro-geodesies. However, we 
show that, in general, it is not possible to take the two counterrotating fluids as circulating along 
electro-geodesies nor take the two counterrotating tangential velocities as equal and opposite. A 
simple family of models of counterrotating charged disks based on the Kerr-Newman solution are 
considered where we obtain some disks with a CRM well behaved. We also show that the disks 
constructed from the Kerr-Newman solution can be interpreted, for all the values of parameters, 
as a matter distribution with currents and purely azimuthal pressure without heat flow. The 
models are constructed using the well-known "displace, cut and reflect" method extended to 
solutions of vacuum Einstein-Maxwell equations. We obtain, in all the cases, counterrotating 
Kerr-Newman disks that are in agreement with all the energy conditions. 
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1 Introduction 



Several methods are known to exactly solve the Einstein and Einstein-Maxwell field equations, or 
to generate new exact solutions from simple known solutions pQ. However, the above mentioned 
methods in general lead to solutions without a clear physical interpretation or to solutions that 
depend on many parameters without a clear physical meaning. Accordingly, it is of importance to 
have some appropriate procedures to obtain physical interpretations of these exact solutions. So, 
in the past years such procedures have been developed for static and stationary axially symmetric 
solutions in terms of thin and, more recently, thick disk models. 

Stationary or static axially symmetric exact solutions of Einstein equations describing relativistic 
thin disks are of great astrophysical importance since they can be used as models of certain stars, 
galaxies and accretion disks. These were first studied by Bonnor and Sackfield |2J, obtaining pres- 
sureless static disks, and then by Morgan and Morgan, obtaining static disks with and without radial 
pressure [31 H]. In connection with gravitational collapse, disks were first studied by Chamorro, Gre- 
gory, and Stewart. Also thin disks with radial tension were considered pa]. Several classes of exact 
solutions of the Einstein field equations corresponding to static thin disks with or without radial 
pressure have been obtained by different authors [7J - [15] . 

Rotating thin disks that can be considered as a source of a Kerr metric were presented by Bicak 
and Ledvinka [IB], while rotating disks with heat flow were studied by Gonzalez and Letelier |17j . 
The nonlinear superposition of a disk and a black hole was first obtained by Lemos and Letelier |lUj . 
Perfect fluid disks with halos were studied by Vogt and Letelier [18J. The stability of some general 
relativistic thin disks models using a first order perturbation of the energy-momentum tensor was 
investigated by Ujevic and Letelier [19]. 

Gonzalez and Letelier [20] constructed models of static relativistic thick disks in various coordinate 
systems. Although the disks have constant thickness, the matter density decreases rapidly with radius 
and the z coordinate, and in principle they also can be used to represent both the disk part and the 
central bulges of galaxies. Also Vogt and Letelier [21] considered more realistic three-dimensional 
models for the gravitational field of Galaxies in the General Relativistic context. Essentially they 
formulate the General Relativistic versions in isotropic coordinates of the potential-density pairs 
deduced by Miyamoto and Nagai [221 [23] and Satoh [2J. 

Disk sources for stationary axially symmetric spacetimes with magnetic fields are also of as- 
trophysical importance mainly in the study of neutron stars, white dwarfs and galaxy formation. 
Although disks with electric fields do not have clear astrophysical importance, their study may be of 
interest in the context of exact solutions. Thin disks have been discussed as sources for Kerr- Newman 
fields [23], magnetostatic axisymmetric fields [26J, conformastationary metrics [27], while models of 
electrovacuum static counterrotating dust disks were presented in [23]. Charged perfect fluid disks 
were also studied by Vogt and Letelier [29], and charged perfect fluid disks as sources of static and 
Taub-NUT-type spacetimes by Garcia- Reyes and Gonzalez [30], [31] • 

In all the above cases, the disks are obtained by an "inverse problem" approach, called by Synge 
the "g-method 1 |32j. The method works as follows: a solution of the vacuum Einstein equations is 
taken, such that there is a discontinuity in the derivatives of the metric tensor on the plane of the disk, 
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and the energy-momentum tensor is obtained from the Einstein equations. The physical properties 
of the matter distribution are then studied by an analysis of the surface energy-momentum tensor 
so obtained. Another approach to generate disks is by solving the Einstein equations given a source 
(energy-momentum tensor). Essentially, they are obtained by solving a Riemann-Hilbert problem 
and are highly nontrivial [33] -[39]. A review of this kind of disks solutions to the Einstein-Maxwell 
equations was presented by Klein in [ID]. 

Now, when the inverse problem approach is used for static electrovacuum spacetimes, the energy- 
momentum tensor is diagonal and its analysis is direct and, except for the dust disks, all the obtained 
disks have anisotropic sources with azimuthal stress different from the radial stress. On the other 
hand, when the considered spacetime is stationary, the obtained energy-momentum tensor is non- 
diagonal and the analysis of its physical content is more involved and, in general, the obtained source 
is not only anisotropic but with nonzero heat flow. Due to this fact, there are very few works about 
of stationary electrovacuum disks and they are limited to disks obtained with solutions that lead to 
disks without heat flow [2S1 EI] • 

The necessary condition to obtain a thin disk without heat flow is that the eigenvalues of the 
energy-momentum tensor must be real quantities, which can be only for very few known electrovac- 
uum solutions. In [3 lj we consider a Taub-NUT type solution such that the energy-momentum 
tensor can be written as an upper right triangular matrix, so that the diagonalization is trivial and 
the eigenvalues are real quantities. However, the obtained disks are not really rotating disks since the 
spatial components of their velocity vectors are zero with respect to the coordinates and so the disks 
are "locally statics". The first true rotating electrovacuum thin disks were obtained by Ledvinka, 
Bicak, and Zofka [25] by applying the "displace, cut and reflect" method to the Kerr-Newman solu- 
tion. The so obtained disks have no radial pressure and no heat flow. However, the authors do not 
show if the eigenvalues of the momentum-energy tensor are real quantities for all the values of the 
parameters, is that, if these disks can always be interpreted as a matter distribution with currents 
and purely azimuthal pressure, or if there are some case where can exist nonzero heat flow (complex 
eigenvalues) . 

The above disks can also be interpreted as made of two counterrotating streams of moving charged 
particles, as was also indicated in [25J. Now, in order to do this interpretation, the counterrotating 
tangential velocities of the two streams must to satisfy a constraint, which in general is not satisfied 
for disks obtained from generic stationary electrovacuum solutions. In [22] the authors take the two 
counterrotating streams as circulating along electro-geodesies, but they do not show if such decom- 
position can be done. In addition, as we will show in this paper, in general the electro-geodesies 
motion do not agree with the above mentioned constraint and so it is necessary to consider another 
possibility for the complete determination of the counterrotating tangential velocities. Another pos- 
sibility, also commonly assumed, is to take the two counterrotating velocities as equal and opposites 
but, as we will show, in general the counterrotating velocities are not completely determined by the 
constraint, so that the corresponding interpretation as two counterrotating streams is not possible. 

The above interpretation is obtained by means of the Counterrotating Model (CRM) in which 
the energy-momentum tensor of the source is expressed as the superposition of two counterrotating 
perfect fluids. Now, even though this interpretation can be seen as merely theoretical, there are 
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observational evidence of disks made of streams of rotating and counterrotating matter (see, for 
instance, [HI H21 S3]). These disks are made of stars and gas so that they are disks with pressure. 
Nevertheless, as is suggested in [42J , the preexisting galaxies have a component originally constituted 
mainly by a gas free stellar disk, i.e., collisionless matter or dust. A detailed study of the CRM for 
generic relativistic static thin disks was presented in [15] for the vacuum case, whereas the extension 
for static electrovacuum disks was presented in [281 [30]. O n the other hand, the CRM for stationary 
thin disks has not been completely developed, neither for the vacuum case, and only a preliminary 
version of it was presented in [17] for the case of stationary thin disks without heat flow and with 
positive radial stress (pressure). 

The purpose of the present paper is twofold. In first instance, we present a detailed analysis of 
the energy-momentum tensor and the surface current density for electrovacuum stationary axially 
symmetric relativistic thin disks of infinite extension without radial stress, in the case when the 
energy-momentum tensor of the disks can be diagonalized, so that there is not heat flow. And, in 
the second place, we present the complete study of the Counterrotating Model for these stationary 
thin disks. The paper is structured as follows. In Sec. II we present a summary of the procedure 
to obtain models of rotating thin disks with a purely azimuthal pressure and currents, using the 
well-known "displace, cut and reflect" method extended to solutions of Einstein-Maxwell equations, 
in the case when the eigenvalues of the energy-momentum tensor of the disk are real quantities. In 
particular, we obtain expressions for the surface energy-momentum tensor and the surface current 
density of the disks. 

In Sec. Ill the disks are interpreted in terms of the counterrotating model (CRM). We find the 
general constraint over the counterrotating tangential velocities needed to cast the surface energy- 
momentum tensor of the disk as the superposition of two counterrotating charged dust fluids. We then 
show that this constraint can be satisfied if we take the two counterrotating tangential velocities as 
equal and opposite as well as by taking the two counterrotating streams as circulating along electro- 
geodesies. However, we show that, in general, it is not possible to take the two counterrotating 
fluids as circulating along electro-geodesies nor take the two counterrotating tangential velocities as 
equal and opposite. We also find explicit expressions for the energy densities, current densities and 
velocities of the two counterrotating fluids. 

In the following section, Sec. IV, we consider a family of models of counterrotating charged dust 
disks based on the Kerr-Newman metric, perhaps the only simple electrovacuum solution that lead 
to stationary thin disks without heat flow. We show that for Kerr-Newman fields the eigenvalues 
of the energy-momentum tensor of the disks are always real quantities, for all the values of the 
parameters, and so they do not present heat flow in any case. We also analyze the CRM for these 
disks and study the tangential velocities, energy and electric charge densities of both streams when 
the two fluids move along electrogeodesics and when they move with equal and opposite velocities. 
Also the stability against radial perturbation is analyzed in both of the cases. Finally, in Sec. V, we 
summarize our main results. 
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2 Electrovacuum rotating relativistic thin disks 

A sufficiently general metric for our purposes can be written as the Weyl-Lewis-Papapetrou line 
clement |1 j 

ds 2 = - e™(dt + Wdipf + e~ 2 VV + e 2A (dr 2 + dz% (1) 

where we use for the coordinates the notation x 1 , x 2 , x 3 ) = (t,<p,r,z), and W, and A are 
functions of r and z only. The vacuum Einstein-Maxwell equations, in geometric units in which 
87rG = c = /io = eo = 1, are given by 

Gab = T a b, (2a) 

F ab ]b = 0, (2b) 

with 

Tab = F ac F b c — -^g ab F cd F cd (3a) 

F a b = A bA — Aa t b, (3b) 

where A a = (A t , A v , 0, 0) and the electromagnetic potentials A t and A v are also functions of r and 
z only 

For the metric (CQ), the Einstein-Maxwell equations are equivalent to the system [44] 



V • [r- 2 f(yA v - W VA t ] = 0, (4a) 

V ■ VA t + r~ 2 fW{VA^ - W VA t ] = 0, (4b) 

V • [r~ 2 f 2 VW - 2r' 2 fA t {WA v - W VA t )} = 0, (4c) 
/V 2 / = V/ • V/ - r- 2 / 4 VW • VW + fVA t ■ VA t 

2 f 3i 



+ r~ 2 f 3 (VA v - WW A t ) ■ (VA^ - WVA t ), (4d) 
i(W 2 -W 2 )e--l 



A, r = r(^ r -^ z )-^(Wl-Wl)e^-^(r 2 e^-W 2 e^)(A 2 -A 2 



+ i(Aj >r - <Je 2 * - £w(iVV - 4**A*)e 2 *, (4e) 
A, z = 2rvl/ ir vl> iZ -lw. r W. 2 e 4 *-i(rV 2 ^-W 2 e 2vI ')A, r A. 2 

+ ^V^e 2 * - ^W(i,, r At, 2 + A^A^e™ , (4f) 

where V is the standard differential operator in cylindrical coordinates and / = e 2 *. 

In order to obtain a solution of (I2~a|) - (j2bjl representing a thin disk at z = 0, we assume that the 
components of the metric tensor are continuous across the disk, but with first derivatives discontin- 
uous on the plane z = 0, with discontinuity functions 

bab = 9ab,z\ z=Q+ ~ Qab,z\ z=Q - = % 9ab,z\ z=0+ - (5) 
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Thus, by using the distributional approach [I6j HH 08] or the junction conditions on the extrinsic 
curvature of thin shells EOl [51] , the Einstein-Maxwell equations yield an energy- momentum tensor 
T ab = T* m + TX\ where T^ at = Qab 5(z), and a current density J a = jj(z) = -2e 2 (*~ A ) A a<z 6(z), 
where S(z) is the usual Dirac function with support on the disk. T^ m is the electromagnetic tensor 
defined in Eq. ( I3al) . j a is the current density on the plane z — 0, and 

Ql = \{b az 8 z h - b zz 5 a b + g a % - g zz b a b + ¥ c (g zz 5 a b - g az 5 z b )} 

is the distributional energy-momentum tensor. The "true" surface energy-momentum tensor (SEMT) 
of the disk, S ab , and the "true" surface current density, j a , can be obtained through the relations 



S ab = J TT ds n = e A ~* Q ab , (6a) 
j a = I Ja ds n = e A ~* j a , (6b) 



where ds n = y/g zz dz is the "physical measure" of length in the direction normal to the disk. 
For the metric ([!]), the nonzero components of S b a are 



Si = ^— [2r 2 (A„ - 2*,, ) - e m WW, z } , (7a) 

S* = -^[4r 2 WV, z + (r 2 + WV>, 2 ], (7b) 

So 1 = ^[e 4 *W„], (7c) 

Si =^[2^+6^], (7d) 



e 

r 

and the nonzero components of the surface current density j a are 

Jt = -2e*- A ^, 2 , (8a) 
j v = -2e*-%, 2 , (8b) 

where all the quantities are evaluated at z = + . 

These disks are essentially of infinite extension. Finite disks can be obtained introducing oblate 
spheroidal coordinates, which are naturally adapted to a disk source, and imposing appropriate 
boundary conditions. These solutions, in the vacuum and static case, correspond to the Morgan 
and Morgan solutions [3]. A more general class of solutions representating finite thin disks can be 
constructed using a method based on the use of conformal transformations and solving a boundary- 
value problem El H51 EB ED] • 

Now, in order to analyze the matter content of the disks is necessary to compute the eigenvalues 
and eigenvectors of the energy-momentum tensor. The eigenvalue problem for the SEMT fl7al) - (17dl) 



s a b e = a r, (9) 
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has the solutions 

A, 



\(t±Vd), (io) 



where 

T = S° + Si, D = (Si - S° ) 2 + 4 S? Si, (11) 
and X r = X z = 0. For the metric (TTT) 



/ = 4e»- A (A,,-*. ; ). (121.) 



D = 4^-^-(4r 2 q 2 z -Wy*) = A 2 -B 2 , (12a) 



where 

A = 4^ z e*" A , S = -W z e 3 *- A . (13) 

r 

The corresponding eigenvectors are 

a = ( e±, ei, o, o), 

X a = e C/-A( , 0,1,0), (14) 
ya = e ^-A(0, 0,0,1), 

with 

^ ) = 2N ^( ss-S±Vd )'' (15) 

where 

N ± = VD(-VD±A). (16) 

We only consider the case when D > 0, so that the two eigenvalues A± are real and different 
and the two eigenvectors are orthogonal, in such a way that one of them is timelike and the other is 
spacelike. Since \A\ > \^D, from (Tl6l) follows that when A > the negative sign corresponds to the 
timelike eigenvector while the positive sign to the spacelike eigenvector. When A < we have the 
opposite case. So the function $ }Z determines the sign of the norm. 

Let V a be the timelike eigenvector, V a V a = —1, and W a the spacelike eigenvector, W a W a = 1. 
In terms of the orthonormal tetrad or comoving observer ea 6 = {V b , W b , X b , Y b }, the SEMT and the 
surface electric current density may be decomposed as 

S ab = eV a V b +p v W a W b , (17a) 

L = fV a +fW a , (17b) 

where 

e = -A±, V v = A T , (18) 
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are, respectively, the surface energy density, the azimuthal pressure, and 



r = -v a j a , f = w a j a , (19) 

are the surface electric charge density and the azimuthal current density of the disk measured by 
this observer. In (1181) the sign is chosen according to which is the timelike eigenvector and which is 
the spacelike eigenvector. However, in order to satisfy the strong energy condition g = e + p v > 0, 
where g is the effective Newtonian density, we must choose as the timelike eigenvector and £ + 
as the spacelike eigenvector. These condition characterizes a disk made of matter with the usual 
gravitational attractive property. Consequently z must be taken positive. So we have 

e = -A_, Pip = X+, (20) 

and 

V° = -^={Sl-S{-VD), (21a) 

= 2ue ~* ci (21b) 

where v = ± so that the sign is chosen according to the causal character of the timelike eigenvector 
(observer's four- velocity), 

W° = -7^=S?, (22a) 
V2M 

W 1 = —^(Sl- S° + y/D), (22b) 
V2M 



where 

M 



Vd{ 9ii Vd + 2rWB + (r 2 e-™ + W 2 e 2 *)A} . (23) 



3 Counterrotating charged dust disks 

We now consider, based on Refs. [15] and [30], the possibility that the SEMT S ab and the current 
density j a can be written as the superposition of two counterrotating charged fluids that circulate in 
opposite directions; that is, we assume 

gab = gab + gab ? 
f = f + +f-: (24b) 

where the quantities on the right-hand side are, respectively, the SEMT and the current density of 
the prograde and retrograde counterrotating fluids. 
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Let U± = (U±, U±, 0, 0) = Z7° (1, u±, 0, 0) be the velocity vectors of the two counterrotating fluids, 
where u± = U±/U± are the angular velocities of each stream. In order to do the decomposition (I24al) 
and (I24bl) we project the velocity vectors onto the tetrad ea b , using the relations [52] 

Ul = e\U b ± , Ul = % a U\. (25) 
In terms of the tetrad (fl4l) we can write 

V a + v+W a 
V 1 - v ± 



so that 



v = VIE^uiWiES^l, (27a) 



^ = yiZS^yTZ^ (27b) 

f + — v_ 

where v± = U±/U± are the tangential velocities of the streams with respect to the tetrad. 

Another quantity related with the counterrotating motion is the specific angular momentum of 
a particle rotating at a radius r, defined as h± = g^JJ%. This quantity can be used to analyze the 
stability of circular orbits of test particles against radial perturbations. The condition of stability, 

d(h 2 ) . , 

is an extension of Rayleigh criteria of stability of a fluid in rest in a gravitational field |53j. For an 
analysis of the stability of a rotating fluid taking into account the collective behavior of the particles 
see for example Refs. [Ml EH] . 

Substituting (I27ap and (127bp in (jl7ajl we obtain 



gab = 
+ 



F(v-,v-)(l-v 2 + ) U%U\ 

{v + — f_) 2 
F(v+,v+)(l -v 2 _) XJ a JJ h _ 

(v + — f_) 2 

F(v+,v-)(l -vl)*(l -v 2 _)^{U%U b _ + UtU b + ) 
(v + — f_) 2 



where 



F(v 1 ,v 2 ) = eu!V 2 +p v . (29) 

Clearly, in order to cast the SEMT in the form (|24aj) . the mixed term must be absent and therefore 
the counterrotating tangential velocities must satisfy the following constraint 

F{y + ,v-) = ev + v- + p v = 0, (30) 
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where we assume that |t>±| ^ 1. 

Then, assuming a given choice for the tangential velocities in agreement with the above relation 
we can write the SEMT as (I24al) with 



(31) 



so that we have two counterrotating dust fluids with surface energy densities, measured in the 
coordinates frames, given by 

r 1 „ ; 2 



e±= ±- ev T , (32) 

y T - V±\ 

Thus the SEMT S ab can be written as the superposition of two counterrotating dust fluids if, and 
only if, the constraint ( |30l) admits a solution such that v + ^ t>_. 

Similarly, substituting ( I27aj) and ( I27bl) in ( 117bl) we can write the current density as ( I24bl) with 



where a± are the surface electric charge densities, measured in the coordinates frames, 



v± - v T 



if 



(33) 



(34) 



Thus, we have a disk makes of two counterrotating charged dust fluids with surface energy densities 
given by fl32l) . and surface electric charge densities given by fl34l . 

As we can see from Eqs. fl26|) . fl32l) and fl34l . all the main physical quantities associated with the 
CRM depend on the counterrotating tangential velocities v±. However, the constraint ( |30l) does not 
determine v± uniquely so that we need to impose some additional requirement in order to obtain a 
complete determination of the tangential velocities leading to a well defined CRM. 

A possibility, commonly assumed [251 [38], is to take the two counterrotating streams as circulat- 
ing along electrogeodesics. Now, if the electrogeodesic equation admits solutions corresponding to 
circular orbits, we can write this equation as 

le± 9ab , r U£Ul = -a±F ra U a ± . (35) 



In terms of u± we obtain 



2 e ±( C/ ±) 2 (fl , ll,r w ± + 2 901,rV± + #00/, 



-a±Ul{A t , r + A Vtr u±). 



From flmj) , fl24bl . (|3T|), and ([33]) we have 



e±iU° ± ) 2 



UJ± — LU T 



LU T S°° 



w 3 



S 11 - m±S 01 
S 01 -u±S 00 ' 



(36) 

(37a) 
(37b) 
(37c) 
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and substituting (137aj) and (137b j) in (136]) we find 



-(S 01 - u; T S ,00 )(5fii ir ^± + 2^oi,rW± + #ocv, 



and using (I37cl) we obtain 

^S 01 ) 2 -S™S u }(gn,rul+2goi,rUJ ± +goo t ', 
Therefore we conclude that 



-(j 1 - <*ff){At, r + A^ r u ± ) 



(3J 



[^Y-^^W-^lCAr+^rWi). (39) 



-iOI ;0 oOO -1 



-T 2 ±JW^Tyr 3 



with 



7i 

r 3 



001,7- + A, 



90Q,r + 2At >r 



I 7 ! 

fS 01 -j l S 00 
v,r S 01 S 01 — S 00 S n ' 

jOgOl.jlgOO 
gOl gOl _ gOOgll 

f S oi_f S n 



fS m -fS n 

gOlgOl _ gOOgll ' 



.4 



' gOlgOl _ gOOgll ' 

On the other hand, in terms of u± we get 

~W + W l u± 

v± = — 

V + V 1 u ± 



T X W% - 2T 2 W W 1 + T 3 W{ 



and so, by using (HO]) , we have that 



V + V- 



2 ' 



so that, using (I17al) . we get 



F(v + ,v. 



32e 4(*-A) A 2 l ( r 2 A ^ + 4r 2 A ^ *-A _ W 2 5*- 



r 3 (A + VD)p^S 01 S m - SWS^frVf - 2T 2 V Vi + T 3 V { 

x A 2 -2r$ r * 2 + -W r W z e w + i(rV 2 * - W 2 e 2 *)A tr ^ 
2r r 

~~A iPjr A iP:Z e + —W(A,p jr At tZ + A^^ z At^)e 



t,z 



2* 



(40) 

(41a) 
(41b) 
(41c) 

(42) 



(43) 



(44) 



Finally, using the Einstein-Maxwell equation ( I4fj) follows immediately that -F(t> + ,t>_) vanishes and 
therefore the electrogeodesic velocities satisfy the constraint (130]) and so, if the electrogeodesic equa- 
tion admits solutions corresponding to circular orbits, we have a well defined CRM. 
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Another possibility is to take the two counterrotating fluids not circulating along electrogeodesics 
but with equal and opposite tangential velocities, 



v ± = ±v = ±y/p tp /e. (45) 

This choice, that imply the existence of additional interactions between the two streams (e.g. col- 
lisions), leads to a complete determination of the velocity vectors. However, this can be made only 
when < \p v /e\ < 1. In the general case, the two counterrotating streams circulate with different 
velocities and we can write (I5U|) as 

v+v- = (46) 
e 

However, this relation does not determine completely the tangential velocities, and therefore the 
CRM is undetermined. 

In summary, the counterrotating tangential velocities can be explicitly determined only if we 
assume some additional relationship between them, like the equal and opposite condition or the 
electro-geodesic condition. Now, can happen that the obtained solutions do not satisfy any of these 
two conditions. That is, the counterrotating velocities are, in general, not completely determined 
by the constraint (1301) . Thus, the CRM is in general undetermined since the counterrotaing energy 
densities and pressures can not be explicitly written without a knowledge of the counterrotating 
tangential velocities. 



4 Disks from a Kerr-Newman solution 

As an example of the above presented formalism, we consider the thin disk models obtained by means 
of the "displace, cut and reflect" method applied to the well known Kerr-Newman solution, which 
can be written as 



where a 2 + b 2 = c 2 > 1, with 
k 



A 
W 

A t 

A, n 



I In 

2 



In 



a 2 x 2 + b 2 y 2 — c 
(ax + c 2 ) 2 + b 2 y 2 



2„,2 



ax + by 



a 2 (x 2 — y 2 ) 



c 2 kb(l - y 2 )(2ax + 1 + c 2 ) 

a(a 2 x 2 + b 2 y 2 — c 2 ) 
cy/2(c 2 - l)(ax + c 2 ) 
(ax + c 2 ) 2 + b 2 y 2 

-k-(l-y 2 )A t , 
a 



1 



m(l — qq*) 



m(l — qq*) 



V 1 - qq*' 



k = Vm 2 — L 2 — 



e 

^ m' 



(47a) 
(47b) 
(47c) 

(47d) 
(47e) 

(48) 
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where m, L and e are the mass, angular momentum and electric charge parameters of the Kerr- 
Newman black hole, respectively. Note that c is the parameter that controls the electromagnetic 
field. The prolate spheroidal coordinates, x and y, are related with the Weyl coordinates by 

r 2 = k 2 {x 2 - 1)(1 -y 2 ), z + z = kxy, (49) 

where 1 < x < oo, < y < 1, and k is an arbitrary constant. Note that we have displaced the 
origin of the z axis in zq. This solution can be generated, in these coordinates, using the well-known 
complex potential formalism proposed by Ernst [B] from the Kerr vacuum solution [1] . When c = 1 
this solution reduces to the Kerr vacuum solution. 
Let be D = k 2 D, T = kT, and j t = kj t , therefore 

~ Ac 2 ay{2x(l - y 2 )(x 2 + lax + c 2 ) - (x 2 - y 2 )[a(x 2 + 1) + x(l + c 2 )]} 

(x 2 -y 2 ) 3 / 2 [(ax + c 2 ) 2 + b 2 y 2 } 3 / 2 ' ( ^ 

D = 16c 4 a 2 y 2 {[a(x 2 + 1) + x(l + c 2 )] 2 - b 2 f 2 } 
(x 2 — y 2 )[(ax + c 2 ) 2 + b 2 y 2 ] 3 



Jt 



2ca/2(c 2 - l)ay{-b 2 y 2 {3ax 2 + 2xc 2 - a) + (ax + c 2 ){a 2 x 3 + ac 2 x 2 

(x 2 - y 2 yi 2 [(ax + c 2 ) 2 + b 2 y 2 f' 2 
-a 2 x + 2b 2 x - ac 2 )}, (50c) 
2cy/2(c 2 - l)by(l - y 2 ){-ab 2 y 2 {x 2 - 1) + (ax + c 2 )(3a 2 x 3 
(x 2 - y 2 y/ 2 [(ax + c 2 ) 2 + b 2 y 2 f/ 2 
+5ac 2 x 2 - a 2 x + 2b 2 x + 2c A x - ac 2 )}. (50d) 



In the above expressions x and y are given by 



2x = ^/r 2 + {a + l) 2 + ^r 2 + [a - l) 2 , (51a) 
2y = y/f 2 + (a + l) 2 - ^/r 2 + (a - l) 2 , (51b) 

where r = r/k and a = Zo/& ; with a > 1. 

Now, in order to analyze the behavior of D, is enough to consider the expression 

D = [a(x 2 + l)+x(l + c 2 )} 2 -b 2 f 2 , (52) 

that can be written as 

D = a(l + c 2 )R + [a(a-l) + 2 + f} + a(l+c 2 )R4a(a + l) + 2 + f] 

+ \R + R.[(c 2 + l) 2 + a 2 (f 2 + a 2 + 3)] + |[(c 2 + l) 2 + 5a 2 ] (53) 

+ lf2 [c 4 + ! + a 2 ( ~2 + 2a 2 + 6)] + l a 2 [(c 2 + 1)2 + a 2 (a2 + 3^ 
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where R± = \/f 2 + (a ± l) 2 . Since a(a =F 1) + 2 > for any a, from (|53|) follows that L> always is 
a positive quantity for Kerr-Newman fields and therefore the eigenvalues of the energy-momentum 
tensor are always real quantities. So we conclude that these disks can be interpreted, for all the values 
of parameters, as a matter distribution with currents and purely azimuthal pressure and without heat 
flow. 

We can see also that for the vacuum case, when c = 1, D is everywhere positive. These disks, 
obtained from the Kerr vacum solution, were previously considered by Gonzalez and Letelier in the 
reference [17]. In this previous work, due to a mistake in the computation of the expressions for T 
and D, was concluded that the energy- momentum tensor could present complex eigenvalues for some 
values of the parameters. As we can see from the expressions presented here, this is not correct and, 
in all the cases, we have a matter distributions with purely azimuthal pressure and without heat flow 
for all the values of parameters. 

In order to study the behavior of the main physical quantities associated with the disks, we 
perform a graphical analysis of them for disks with a = 2, b = 0.2 and c = 1.0, 1.5, 2.0, 2.5, and 
3.5, as functions of f. For these values of the parameters we find that ^ z is a positive quantity in 
agreement with the strong energy condition. Therefore, e = — A_ and = X + . However, one also 
finds values of the parameter for which ty z takes negative values. Furthermore, we take v — — 1 in 
order to V a be a future-oriented timelike vector. 

In Fig. [T]we show the surface energy density e and the azimuthal pressure p^. We see that the 
energy density presents a maximum at f = and then decreases rapidly with f, being always a 
positive quantity in agreement with the weak energy condition. We also see that the presence of 
electromagnetic field decreases the energy density at the central region of the disk and later increases 
it. We can observe that the pressure increases rapidly as one moves away from the disk center, reaches 
a maximum and later decreases rapidly. We also observe that the electromagnetic field decreases the 
pressure everywhere on the disk. 

The electric charge density j t and the azimuthal current density j , measured in the coordinates 

frame, are represented in Fig. [2} whereas the electric charge density j and the azimuthal current 

density j , measured by the comoving observer, are represented in Fig. |3j We observe that the 
electric charge density has a similar behavior to the energy and that the current density have a 
similar behavior to the pressure which is consistent with the fact that the mass is more concentrated 
in the disks center. We also computed this functions for other values of the parameters and, in all 
the cases, we found the same behavior. 

We now consider the CRM for the same values of the parameters. We first consider the two 
counterrotating streams circulating along electrogeodesics. In Fig. H] we plot the tangential velocity 
curves, v + and t>_. We see that these velocities are always less than the light velocity. We also see 
that the inclusion of the electromagnetic field make less relativistic these disks. In Fig. [5] we have 
drawn the specific angular momenta h 2 + and h 2 _ for the same values of the parameters. We see that 
the presence of electromagnetic field can make unstable these orbits against radial perturbations. 
Thus the CRM cannot apply for c = 6 (bottom curve). In Fig. Owe have plotted the surface energy 
densities e + and e_. We see that these quantities have a similar behavior to the energy density e. 
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In Fig. [7J we plotted the surface electric charge densities a + and <r_. We find that these quantities 
have also a similar behavior to e±. 

The Figs. H]- [7] show that the two fluids are continuous in r which implies to have two particles in 
counterrotating movement in the same point in spacetime. So this model could be possible when the 
distance between streams (or between the counterrotating particles) were very small in comparing 
with the length r so that we can consider, in principle, the fluids continuous like is the case of 
counterrotating gas disks present in disk galaxies. 

Finally, in the case when the two fluids move with equal and opposite tangential velocities (non- 
electrogeodesic motion) we find that the physical quantities have a similar behavior to the previous 
one. 

5 Discussion 

We presented a detailed analysis of the energy-momentum tensor and the surface current density 
for electrovacuum stationary axially symmetric relativistic thin disks of infinite extension without 
radial stress, in the case when the energy-momentum tensor of the disks can be diagonalized, so that 
there is not heat flow. The surface energy-momentum tensor and the surface current density were 
expressed in terms of the comoving tetrad and explicit expressions were obtained for the kinematical 
and dynamical variables that characterize the disks. That is, we obtained expressions for the velocity 
vector of the disks, as well for the energy density, azimuthal pressure, electric charge density and 
azimuthal current density. 

We also presented in this paper the stationary generalization of the Counterrotating Model (CRM) 
for electrovacuum thin disks previously analyzed for the static case in [281 EOJ. Thus then, we were 
able to obtain explicit expressions for all the quantities involved in the CRM that are fulfilled when 
do not exists heat flow and when we do not have radial pressure. We considered both counter rotation 
with equal and opposite velocities and counter rotation along electrogeodesics and, in both of the 
cases, we found the necessary conditions for the existence of a well defined CRM. 

A general constraint over the counterrotating tangential velocities was obtained, needed to cast the 
surface energy-momentum tensor of the disk in such a way that can be interpreted as the superposition 
of two counterrotating dust fluids. The constraint obtained is the generalization of the obtained 
for the vacuum case in [17], for disks without radial pressure or heat flow, where we only consider 
counterrotating fluids circulating along geodesies. We also found that, in general, there is not possible 
to take the two counterrotating tangential velocities as equal and opposite neither take the two 
counterrotating fluids as circulating along geodesies. 

A simple family of models of counterrotating charged disks based on the Kerr-Newman solution 
were considered where we obtain some disks with a CRM well behaved. We also find that the disks 
constructed from the Kerr-Newman solution can be interpreted, for all the values of parameters, as a 
matter distribution with currents and purely azimuthal pressure and without heat flow. We obtain, 
for all the values of parameters, counterrotating Kerr-Newman disks that are in agreement with 
all the energy conditions. Finally, the generalization of these models to the case of electrovacuum 
stationary axially symmetric solutions where the energy-momentum tensor of the disk can to present 
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complex eigenvalues for some values of the parameters, the stability of counterrotating fluids taking 
into account the collective behavior of the particles, and a thermodynamic analysis of the disks, will 
be considered in future works. 
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Figure 1: (a) The surface energy density e and (b) the azimuthal pressure p v for Kerr- Newman disks 
with a = 2, b = 0.2 and c = 1.0 (curves with crosses), 1.5, 2.0, 2.5, 3.0, and 3.5 (bottom curves), as 
functions of f. 
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(a) (b) 

Figure 2: (a) The surface electric charge density j t and (b) the azimuthal current density j for 
Kerr-Newman disks with a = 2, b = 0.2 and c = 1.0 (axis f), 1.5 (top curves), 2.0, 2.5, 3.0, and 3.5 
(bottom curves), as functions of f. 
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Figure 3: For Kerr-Newman disks we plot, as function of r, (a) j with a = 2, b = 0.2, and c = 1.0 
(axis f), 1.5 (bottom curve, away from the center of disk), 2.0, 2.5, 3.0, and 3.5 (top curve, away 
from the center of disk), and (b) j 1 also with a = 2, b = 0.2, and c = 1.0 (axis f), 1.5 (top curve), 
2.0, 2.5, 3.0, and 3.5 (bottom curve). 




Figure 4: The tangential velocities (a) v + and (b) V- for electrogeodesic Kerr-Newman disks with 
a = 2, b = 0.2 and c = 1.0 (curves with crosses), 1.5, 2.0, 2.5, 3.0, and 3.5 (bottom curves), as 
functions of f. 
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Figure 5: The specific angular momenta (a) h\ and (b) h 2 __ for electrogeodesic Kerr- Newman disks 
with a = 2, b = 0.2 and c = f.O (curves with crosses), 1.5, 2.0, 2.5, 3.0, 3.5, and 6.0 (bottom curves), 
as functions of f. 
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Figure 6: The surface energy densities (a) e + and (6) e_ for electrogeodesic Kerr-Newman disks with 
a = 2, b = 0.2 and c = 1.0 (curves with crosses), 1.5, 2.0, 2.5, 3.0, and 3.5 (bottom curves), as 
functions of f. 
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(a) (6) 

Figure 7: The surface electric charge densities (a) a + and (6) 5"_ for electrogeodesic Kerr-Newman 
disks with a = 2, b = 0.2 and c = 1.0 (axis f), 1.5, 2.0, 2.5, 3.0, and 3.5 (top curves away from the 
center of the disk), as functions of f. 
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